IFIC/95-62 
February 1, 2008 



Searching for Be mesons in the ATLAS experiment at LHC* 

M. Galdon", R. Perez Ochoa''''= , M. A. Sanchis-Lozano'^'''t and J.A. Valls''''^ 

(a) Departamento de Fisica Teorica 

(b) Institute de Fi'sica Corpuscular (IFIC) 

(c) Departamento de Fi'sica Atomica, Molecular y Nuclear 
Dr. Moliner 50, E-46100 Burjassot, Valencia (Spain) 



We discuss the feasibility of the observation of the signal from mesons in the ATLAS 
experiment at the LHC. In particular, we address the decay mode Bc—^J/ipn followed by 
the leptonic decay J/ilj—>-fi'^fi~, which should permit an accurate measurement of the 
mass. We performed a Monte Carlo study of the signal and background concluding that 
a precision of ~ 1 MeV for the B^ mass could be achieved after one year of running at 
"low" luminosity. The semileptonic decay Bc^J /ip fi^i'^ is also considered for a possible 
extraction of \Vcb\- 



1. Introduction 

There is a general consensus in the scientific 
community j^] that the scope of a future high- 
luminosity, high-energy hadron collider like LHC 
should not be restricted to the hunting of the 
standard model Higgs and its extensions, or the 
search for supersymmetry. Other topics requiring 
lower luminosities like top and beauty physics de- 
serve in their own right a close attention. 

In particular, a lot of work has been recently 
devoted to the observation of the Be meson re- 
garding both the hadroproduction and decay 
|3| [Q (see references therein). Specifically, we 
shall focus on the feasibility of its detection in 
the ATLAS experiment at the LHC through some 
decay modes. 
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2. Be signal 

At the center-of-mass energy -^s = 14 TeV, the 
cross-section for beauty production is assumed to 
be 500 /lb leading to 5x10^^ bb pairs per year- 
run (10^ s) at a luminosity of £ w oo/^^ cm^^ 
s"^, corresponding to an integrated luminosity of 
^ 10 fb^"'^. The number of bottom pairs reduces, 
however, to 2.3x10^*^ by requiring events with a 
triggering muon coming from either a 6 or a 5 
under the kinematic cuts p± > 6 GeV/c and |?7| < 
2.2 §. 

On the other hand, assuming that the 6-quark 
fragmentation yields a i?c or a B* with probabil- 
ity of the order of 10"'^ the yield of Be mesons 
(not yet triggered) per year of running would be 
roughly ~ 10^". 

Be ^ J/ip TT channel 

This exclusive channel followed by the leptonic 
decay of the J/ip resonance into a pair of oppo- 
sitely charged muons offers several important ad- 



vantages. First of all, it allows for the mass recon- 
struction of the Be meson. Observe also that any- 
one of the two muons can trigger the decay. Be- 
sides, it is very clean topologically with a common 
secondary vertex for all three charged particles, 
two of them (the muons) with the additional con- 
straint of their invariant mass compatible with a 
J/ip. Furthermore, the expected branching frac- 
tion is not too small, about 0.2% which com- 
bined with the branching fraction of the leptonic 
decay of the resonance BR{J/tp — > fJ.^^^) — 6% 
0, yields an overall branching fraction for the 
signal of 10""'. Thus, the number of such events 
turns out to be ~ 10^ per year of running. 

A study has been performed in order to es- 
timate the signal detection efficiency and back- 
ground for the Be J/-0(— > tt channel. 
The Monte Carlo employed for the signal sim- 
ulation corresponds to a sample of Be — > J/tpT^ 
events while for the background we have used a 
sample of inclusive b muon decays generated in 
all the cases with PYTHIA 5.7 . 

Two types of background were considered: ^ 

a) Combinatorial background due to muons 
from semileptonic decays of bb pairs pro- 
duced at the main interaction. Cascade 
contributions such as b^c^^ arc included 
as well for random combinations with any 
other muon in the same event. 

b) Contamination from prompt J/tJj^s in com- 
bination with another charged hadron (in- 
terpreted as a pion) from the main vertex. 
(In fact data released by Tevatron on the 
J/V' yield point out a production rate quite 
larger than initially expected ^.) Incorrect 
tracking may give rise to the reconstruc- 
tion of a (fake) secondary vertex, becom- 
ing a potential source of a large amount of 
background. 

In a first step, we imposed the following cuts 
on events based on kinematic constraints: 

^Muons coming from semileptonic decays of long-lived 
particles such as pions or kaons contribute in a negligi- 
ble amount to trigger rates. On the other hand B decays 
into J/i/i and a charged particle would give an invariant 
mass quite below the Be mass so they are of no concern 



• P±min{trig.fi) = 6 GeV/c ; 
I Vrnax {trig.fi) 1= 2.2 

• PJ_mm(Ai) = 3 GeV/c ; | rjmaxitJ') 1= 2.5 

• PJ_mm(7r) = 1 GeV/c ; | rimax{T^) 1= 2.5 
. M^+^- = Mj,^ ± 50 MeV 

The two first cuts correspond to the requirement 
of the Ist-level B physics trigger leading in our 
case to an efficiency of ^ 15% in triggering one 
of the two muons from the J/ij]- We next take 
into account the detection efhciency for the sig- 
nal after applying the rest of p±^ and ry cuts which 
turns out to be ~ 21%. Setting the efficiency for 
muon identification as 80% and the track recon- 
struction as 95% dl we get a combined detection 
efficiency of ~ 2% leading to an observable signal 
of about 20, 000 events per year of running. 

The last of the cuts described above constrains 
the two muons invariant mass to be compati- 
ble (within two standard deviations [^) with the 
nominal J /ip mass, thus drastically reducing ran- 
dom combinations. However, background of class 
b) can potentially pass all the kinematic cuts by 
a large amount, so another type of rejection is 
required. 

To this end, we adapted to our needs the ver- 
tex reconstruction (i.e. vertex finding and fit- 
ting) routines of the LEP experiment DELPHI 
at CERN [|l0|. The vertex fitting algorithm pro- 
vides as output the coordinates of the secondary 
vertex, the track momenta re-evaluated with the 
vertex constraint and the goodness of the fit by 
means of the total as well as the contribution of 
each single track to it. In particular, we employed 
for background rejection the three spatial coordi- 
nates and the for each fitted secondary vertex 
formed by the two muons and the charged hadron 
(assumed to be a pion) satisfying the above kine- 
matic constraints. The distance between the re- 
constructed vertex and the primary {pp) interac- 
tion point was thereby determined. We shall re- 
fer to it as the decay length even for background 
events of class b). 

Hence, candidate (either signal or background) 



events were required to pass the following extra 
cuts: 

• total < Xo 

2 

• Xi < ^ for each single track-i 

• decay length larger than Lq 

where Xo^ Lq have to be optimized to remove 
the background as much as possible but with a 
good acceptance for the signal. In our analysis 
we found Xo = 8 and Lq = 350 /^m. 

Figure 1 shows the reconstructed )j/^ tt 

mass distribution for the expected signal above 
the surviving background once all the cuts have 
been applied, for an integrated luminosity of 10 
fb-i |. 

In summary, we have found that the self- 
triggering weak decay Bc—^J/tp tt, followed by 
the leptonic decay of the J/^ into two muons, 
could be clearly observed in the ATLAS detector 
at LHC. Under rather conservative assumptions, 
a total number of sa 10, 000 signal events could 
be fully reconstructed after one year run, corre- 
sponding to 10 fb~^ at "low" luminosity (« 10'^^ 
cm~^ s"^). This represents a signal to back- 
ground ratio of about 0.5 with a statistical signifi- 
cance of w 20 standard deviations above a nearby 
almost flat background. The foreseen mass res- 
olution (standard deviation) of the Be meson is 
about 40 MeV representing a precision of about 
1 MeV for the mass measurement. 

Be ^ J/ip /i^ i^fj. channel 

In spite of the fact that this channel does not 
permit the measurement of the Be mass, its sig- 
nature would be quite clean experimentally when 
the J ftp decays into a pair of muons, providing 
a three muon vertex. In the following we shall 
argue that in LHC experiments this decay mode 
could provide a reliable extraction of the mixing 
matrix element | Vet \ ■ 

The expected semileptonic branching ratio 
~ 2% ^ together with the leptonic branching 
ratio for the J/'0 into two muons ~ 6% 0, yields 
an overall branching fraction of order 10^"^. We 
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Figure 1. Reconstructed (/i"*"/!") j/^ tt mass dis- 
tribution after cuts. The nominal value for the 
Be mass was set equal to 6.3 GeV. The solid line 
corresponds to a linear -|-Gaussian fit 



next combine the last BF with the acceptance of 
the muon trigger ~ 15% (where the three possible 
triggering muons per decay have been taken into 
account.) Besides, the detection of the two re- 
maining particles (i.e. the non-triggering muons 
satisfying some p± and ry cuts ) amounts to an 
acceptance for signal events of ~ 12%. Finally, we 
assume an identification efficiency for each muon 
of 0.8 1^ yielding a combined value of 0.8^ ~ 0.5. 

Thereby, the expected number of useful Be de- 
cays could reach several 10^ per year at low 
luminosity. 

Such large statistics together with the fore- 
seen precision in the muon momentum measure- 
ment should permit the experimental access to 
the kinematic region near zero-recoil of charmo- 
nium (with respect to the Be meson) with a good 
energy/momentum resolution. Moreover, very 
stringent cuts can be put on events and hence, 
expectedly, background could be almost entirely 
removed from the event sample permitting a very 

* Semi-electronic Bc^J/^e'^u decays, followed by the J/ip 
decay into fi'^ fJ.~ can be considered as well. In addition, if 
the muon from the semi-muonic decay (1) gives the trigger, 
the event can be selected with the mode J/i/)— > e+e". The 
trigger acceptance for the former is ~ 12% whereas for the 
latter is ~ 3% ||l 



precise determination of the Be lifetime. We ex- 
plore all these issues in the following section. 

Theory inputs and determination of \Vcb\ 

Next we examine in detail the possibility of 
extracting |Vc&| from experimental data. Let us 
start by recalling that 



dw 

1 dBr{Bc^J/il) fi+v) 

T 



dw 



(1) 



The differential branching ratio dBr/dw can be 
written as 

dBr{Bc~^J/tlj fi+i') 



N{B, 



dw 

dN{B^^J/il} fi+iy) 
dw 



(2) 



where dN/dw represents the number of semilep- 
tonic decays per unit of w for a certain integrated 
luminosity, once corrected by detection efficiency; 
N{Bc) is the corresponding total yield of B^s 
produced for such integrated luminosity. (Notice 
that T could be determined accurately from the 
same collected sample of semileptonic Be decays.) 

In order to get N{Bc) it will be convenient to 
compare the Be production rate to the prompt 
■0' yield in pp collisions, for the same integrated 
luminosity. Therefore 



N{B,) 



ajpp^Bc + X) 
a{pp^ip' + X) 



" e^,■BF{'i|;'^^l+^x-) 

where N (jp' ^ ^i'^ fi^ ) stands for the number of 
prompt Tp' to be experimentally detected through 
the muonic decay mode and e^;,' denotes its detec- 
tion efhciency. Contamination from weak decays 
of B mesons into ■0' states should be efficiently 
removed by means of the vertex capability of the 
inner detector. 

Therefore, we suggest normalizing the B^ sam- 
ple with the aid of the ip' yield through expression 



(3). The V' state is probably preferable to the 
J/ip mainly because the former expectedly should 
not be fed down by higher charmonium states 0. 
Instead, data released by Tevatron shows that the 
J/f/; indirect production through Xc intermediate 
states may be also important. Let us also remark 
that the production of the ip' resonance in pp col- 
lisions is more than one order of magnitude larger 
than initially expected |9). 

At sufficiently large p± Q one reasonably ex- 
pects that the main contribution to heavy quarko- 
nia production comes from the splitting of a 
heavy quark or a gluon. Indeed, even though 
the fragmentation process is of higher order in 
as than the "conventional" leading order dia- 
grams I p^ , the former is enhanced by powers of 
Vi-I^Q relative to the latter |^ Although 
the situation is still controversial in the litera- 
ture, there are strong indications from very de- 
tailed calculations (see Q ||T^ and references 
therein) that fragmentation indeed dominates for 
large enough p±^ . In the following, we first exam- 
ine some color-singlet mechanisms, expected to 
contribute largely to the high-p^ inclusive pro- 
duction of prompt heavy quarkonia. 

On the other hand, according to the set of pa- 
pers in perturbative QCD can provide a reli- 
able calculation for the fragmentation function of 
a high-pj^ parton into Be states with only few in- 
put parameters: the QCD running constant as, 
the charm and bottom masses, and the square 
of the radial wave function at the origin |i?(0)|^. 
Note that the largest uncertainty for normaliza- 
tion purposes comes from the (third power of the) 
charm mass in the fragmentation function (see 
for explicit expressions). However, we are only in- 
terested in the ratio 



aipp^Bc + X) 
aipp^Tp' + X) 



(4) 



""See Ref. |12| for an alternative explanation of the ob- 
served surplus found in Tevatron. However this mech- 
anism would require an uncomfortable large radiative BF 
for higher resonances 

^We find an average px ~ 20 GcV/c for Be mesons simply 
passing the kinematics cuts on px ^'Hd I*? I of the decay 
muons 



where each cross section can be written as a con- 
volution of the parton distribution functions of 
the colHding protons, the cross section for the 
hard subprocess leading to the fragmenting gluon 
or heavy quark and the respective fragmentation 
function [Q. Thus notice that the fragmentation 
functions D-^_^g^, Dg^B^ and Dc^^', Dg^^> [|l8|, 
which themselves are independent of the parton- 
level subprocesses, appear combined as a ratio. 
Therefore, those uncertainties coming from the 
common overall factor automatically cancel 
each other. Moreover, those uncertainties intro- 
duced by the parton distribution functions also 
should significantly diminish in r as well. Still 
one must evaluate the ratio 



\RbA0)\' 



(5) 



In fact there is the possibility of expressing the 
above ratio (aside trivial factors) according to the 
general factorization analysis of |]l9| in a more 
rigorous way as 



<0|Of=(iS'o)|0> 



(6) 



<0|Of (3^1)10 > 

where are local four fcrmion operators^ 

and the matrix elements can be evaluated from 
NRQCD. A similar expression holds for the B*. 
(Moreover, the denominator can be determined 
from the measured leptonic width of the ■0'.) 

Recently, an additional color-octet fragmenta- 
tion mechanism [ pO| has been suggested in order 
to reconcile the experimental results on inclu- 
sive 4'' production with theoretical predictions. 
This mechanism assumes the creation from gluon 
fragmentation of a cc pair in a color-octet state, 
in analogy to Xc production jl^ . Then a new 
nonperturbative parameter ^g^^/j is required, ap- 
pearing in r combined as the dimensionless factor 



'*8(i/j' 



<o\ot'es,)\o> 

<0|Of (35i)|0> 



ii?^'(o)r 



(7) 



'^This implies the application of NRQCD to the evaluation 
of the non-perturvative part of the fragmentation function 
of a 6 quark into a be bound state. Therefore four quark 
operators involve both b- and c-quark fields in this case|25|]. 



where we remark that it can be computed in 
terms of NRQCD matrix elements ||l§ ^ 

Furthermore, if we do not neglect the contri- 
bution to the total Be production of those orbital 
excitations like P-wave states, new fragmentation 
functions -Dj_^^^(.p^ come into play Accord- 
ingly, some new Ki factors will appear which can 
be conveniently expressed, for instance in analogy 
to Eqs. (6) and (7), as: 



1 <0|O^(Pj)|0> 



< 0|Of'=(iS'o)|0 > 



l(6c) 



and 



m 



\RbM\' 



l_ <0|Of-(3^i)|0> 
2 <0|0^(i^o)|0> 



(8) 



oc 



8(6c) 



i^s.(o)r 



(9) 



where rh is the reduced mass of the 6 and c quarks. 
The long-distance Hi parameter is related to 
the square of the derivative of the color-singlet 
wave function at the origin. 

Higher Be resonances like D-wa.ve states 
might be further taken into account, introducing 
new nonperturbative parameters involving higher 
derivatives of the wave function at the origin, or 
the NRQCD matrix element analogues. 

In sum, the evaluation of such Ki parameters 
would allow one to complete the computation of 
the ratio r and thereby to obtain from the exper- 
imental measurement of the ^' yield, the total 
number of Be events for an integrated luminosity 
with the aid of expression (3). 

On the other hand, the hadronic transi- 
tion Bc—^J/4' involves two heavy-heavy systems 
whereby heavy quark spin symmetry should be 
still valid as a first order approach m pS]. This 



^^8(i/)') '^^'^ phenomenologically estimated from B or 



T decays as well [L 
cording to a nonre! 



_ |21| |22|. On the other hand, ac- 
ativistic quark model H'g is related to 



a fictitious color-octet wave function at the origin 



permits the introduction of the analogue of the 
Isgur-Wise function for transitions involving dou- 
bly heavy hadrons, denoted by r?i2(wi-i'2)- At the 
non-recoil point vi = V2, we shall write 

< J/^\ A^' \B, > = 2 7712(1) ^TK^e*/ (10) 

where A'^ = cj^-f^b stands for the axial- vector 
current and £3 represents the four-vector polar- 
ization of the J/ip- 

Following similar steps as in the B decay into 



D* 1 26 1, we write 



lim,. 



1 dr{B^^j/tfj fi-^ 



{ 

47r3 



1 



dw 



m^f m\Tl^{\)\V,^\' (11) 



It is of key importance from the theoretical side 
in our proposal, the existence of a single form fac- 
tor in the above expression to be determined the- 
oretically in a rigorous manner. Note that ?7i2(l) 
may be interpreted in first approximation as the 
overlap of the initial and final hadron wave func- 
tions Q 1^^. Moreover, due to the intrinsic non- 
relativistic nature of its heavy constituents, non- 
relativistic QCD on the lattice or even refined po- 
tential based models of hadrons, could be applied 
in order to get t]yi and its slope at zero recoil in 
a reliable way. 

This fact in conjunction with the presumed ac- 
curacy of the J lt\} momentum measurement sug- 
gests determining H4b| in a similar way as in the 
method proposed by Neubert |^ through the 
B^D*ev decay. 

For the purpose of illustration, we have 
depicted in figure 2 some expected typi- 
cal points and error bars for the decay 
Bc^J l^{-^pL~^ ^i^v for a definite prediction 
of the ratio of production cross sections r in 
Eq.(4). 

Under this condition, let us remark that the 
uncertainty on w (i.e the horizontal error bars) 
is of systematic nature, mainly coming from the 
Be flight-path reconstruction. With regard to the 
vertical coordinate rii2{w)-\Vcb\, we have assumed 
that its uncertainty is essentially statistical p|. 
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Figure 2. Hypothetical rii2{w)-\Vcb\ distributions 
of the decay Bc^J/ipfi^i' for a 500 events sample 
tentatively assigned to an integrated luminosity 
of 20 fb~^ (i.e. two years of running at low lumi- 
nosity). It was assumed an spatial resolution for 
the Be flight-path reconstruction of cr « 100 /xm 
and its lifetime is allowed to vary in the inter- 
val T = 0.5 — 1 ps. Horizontal bars show sys- 
tematic errors basically limiting the number of 
experimental points. Vertical bars take into ac- 
count the statistical fluctuations of the sampled 
events in the bin (distributed according to the ex- 
pected differential rate for a pseudoscalar to vec- 
tor semileptonic decay [p6|). Points line up with 
ri'12 set equal to —2. For more details see Ref. ||2^ 



It is important to stress that the quantity to be 
experimentally measured with a large accuracy 
from Bc-^J /tj} jX^v decays can be written as the 
linear combination: 

X 7712(1) Xl^efcl 

In summary, we suggest to keep an open mind 
on the possibility of an alternative determination 
of \Vcb\ from semileptonic decays of Be mesons at 
LHC. Even if at the time LHC will start to run, 
B factories would have already provided a (still) 
more precise value of \Vcb\ than at present via the 
semileptonic B decay, an independent determina- 
tion of it should bring a valuable cross-check. On 
the other hand, a lot of activity is being devoted 
to the analysis of inclusive production of prompt 
charmonium resonances at Fermilab. Therefore 
we conclude that if the fragmentation approach 
to describe the production of heavy quarkonia at 
the large transverse momentum domain is accu- 
rately tested (and tuned) from experimental data. 
Be semileptonic decays could be competitive with 
the B ones. 

Alternatively, one can turn the question round 
and consider \ Vcb\ as a well-known parameter thus 
verifying QCD calculations, either at the produc- 
tion or at the weak decay level. Indeed, frag- 
mentation into heavy quarkonia offers an inter- 
esting check of perturbative QCD and a deep in- 
sight into the nonperturbative dynamics in the 
hadronic formation. Besides, a precise knowledge 
of the Be production rate is a necessary condition 
for the experimental measurement of the absolute 
branching fraction of any of its decay modes. 

REFERENCES 

1. G. Carboni et al. Proceedings of the ECFA 
Large Hadron Collider Workshop, Aachen 
(1990) 

2. K. Cheung and T.C. Yuan, Phys. Lett. B325 
(1994) 481 ; E. Braaten, K. Cheung and T.C. 
Yuan, Phys. Rev. D48 (1993) 5049; D48 
(1993) R5049 ; K. Cheung, Phys. Rev. Lett. 
71 (1993) 3413 

3. F. Albiol et al. ATLAS internal note 
PHYS-NO-0 58 (1994); IFIC/95-24 
ph/9506306| ) 



10. 



11, 



12, 



13 



M.A. Sanchis-Lozano, Nuc. Phys. B440 
(1995) 251 

ATLAS Technical Proposal, CERN/LHCC 
94-43 

M. Lusignoh and M. Masetti, Z. Phys. C51 
(1991) 549 

Particle Data Group, Phys. Rev. D50 (1994) 
M.A. Sanchis-Lozano et al, ATLAS internal 
note, in preparation 

M. Cacciari and M. Greco, Phys. Rev. Lett. 
73 (1994) 1586 ; E. Braaten, M.A. Doncheski, 
S. Fleming and M.L. Mangano, Phys. Lett. 
B333 (1994) 548 

P. Billoir and S. Qian, Nucl. Instr. and Meth. 
A311 (1992) 139 

C-H Chang and Y-Q Chen, Phys. Rev. D49 
(1994) 3399 

Phys. Lett. 
K. Sridhar, 



14 
15 

16 

17 



18 



19 



21 
22 



F.E. Close, 
; D.P. Roy, 
(1995) 537 
E. Braaten, 
(1994) 3176 



B342 (1995) 369 
Phys. Lett. B345 



, Phys. Rev. D50 
Phys. Lett. B332 
E. Braaten ct al., 



B339 



T.C. Yuan 
; J.P. Ma, 
(1994) 398 (1994) 1586 ; 
Phys. Lett. B333 (1994) 548 
G.A. Schuler, CERN-TH.7170/94 

D. P. Roy, K. Sridhar, Phys. Lett 
(1994) 141 

K. Kolodziej, A. Leike and R. Riickl, MPI- 
PhT/95-36 

M. Cacciari, M. Greco, M. L. Mangano and A. 
PetreUi, Phy s. Lett. B356 (19 95) 553, CERN- 
TH/95-129 ( |hep-ph/9505"379| ) 

E. Braaten and T.C. Yuan, Phys. Rev. Lett. 
71 (1993) 1673 

G. T. Bodwin, E. Braaten and G.P. Lepage, 
Phys. Rev. D51 (1995) 1125 

20. E. Braaten and S. Fleming, Phys. Rev. Lett. 
74 (1995) 3327 ; P. Cho and A. Leibovich, 
CALT-68-1988 

H. D. Trotticr, Phys. Lett. B320 (1994) 145 
G.T. Bodwin, E. Braaten, T.C. Yuan and 
G.P. Lepage, Phys. Rev. D46 (1992) R3703 
; G.T. Bodwin, E. Braaten and G.P. Lepage, 
Phys. Rev. D46 (1992) R1914 



23. K. Cheung a nd T.C. Yuan, CPP-94-37 (|hep- 



(hcp- 



ph/9502250D 



24. M.A. Sanchis, Phys. Lett. B312 (1993) 333 
Z. Phys. C62 (1994) 271 



25. T. Mannel and G. A. Schuler, CERN preprint 
CERN-TH.7468/94 

26. M. Neubert, Phys. Rep. 245 (1994) 

27. M. Neubert, Phys. Lett. B264 (199 1) 455 ; 
CERN-TH/95-107 (lhep-ph/950523^ ) 

28. M. Galdon a nd M.A. Sanchis -Lozano, 
IFIC/95-31 ( |hcp-ph/9506307D to appear in Z. 
Phys. C 



